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Abstract. - Recent experimental studies on the stepwize motion of biological molecular motors 
have revealed that the "characteristic distance" of a step is usually less than the actual step size. 
This observation implies that the detailed-balance condition for kinetic rates of steps is violated 
in these motors. In this letter, in order to clarify the significance of the characteristic distance, we 
study a Langevin model of a molecular motor with a hidden degree of freedom. We find that the 
ratio of the characteristic distance to the step size is equal to unity if the dominant paths in the 
state space are one dimensional, while it deviates from unity if the dominant paths are branched. 
Therefore, this parameter can be utilized to determine the reversibility of a motor even under a 
restricted observation. 
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' Single-molecule measurement techniques have expanded 
the possibilities for studying biological macromolecules 
from physical points of view. For instance, molecu- 
lar motors, which move along protein filaments in a 
stepwise manner fueled by the hydrolysis of adenosine 
tri-phosphate (ATP), have been investigated extensively. 
Through many studies, several of the important properties 
of molecular motors, including their step sizes and the ki- 
netic rates of their forward and backward steps, have been 
Experimentally determined. 

' In this letter, we study the modulation of the kinetic 
r,ates by an external force. Such phenomena have been 
examined in recent experiments [1-5]. Interestingly, the 
results of these experiments suggest that the detailed- 
balance (DB) condition for the kinetic rates is violated 
ifi the stepwise motion of a motor. 

' Here, we start from the DB condition for a system that 
Ijas several discrete states [6]. Let mi and m2 denote two 
of these states, and let k2i and fci2 represent the transition 
rates from mi to u\2 and vice versa, respectively. Then, 
the DB condition is given as 



^=exp(/3AGi2), 

fcl2 



(1) 



where /3 = {k-QT)~^ is the inverse temperature of the heat 
bath (fee being the Boltzmann constant), and AG12 = 
Gi — G2 represents the difference in the (free) energy be- 
tween the states mi and m2. 

Note that the DB condition given in eq. (1) is de- 
rived from a more microscopic condition termed the lo- 
cal detailed-balance (LDB) condition with several assump- 
tions. Let ri and r2 denote microstates in the states mi 
and m2, respectively, and let [r] represent a microscopic 
path connecting ri and r2 in an interval t. Let P([r]|ri) 



represent the transition probability of a path [r] with given 
ri. Then, the LDB condition is expressed as 

P([r]|ri^ 



P{[r]\r2) 



exp [/3AG(ri,r2)] 



(2) 



where [r] is the time- reversed path of [r], and 
AG(ri,r2) = G(ri) — G(r2) represents the difi'erence 
in the (free) energy between the microstates ri and r2. 
Eq. (1) is derived from eq. (2) provided that the ther- 
mal energy is sufficiently less than the hight of a barrier 
separating mi and m2, and the curvature of G(r) around 
these states are the same. It has been known that the 
LDB condition is an essential condition in order to make 
a stochastic model compatible with thermodynamics, and 
to derive several nonequilibrium equalities, including the 
fluctuation theorems and the Jarzynski equality [7,8]. Fur- 
thermore, it has been argued that eq. (2) is satisfied in a 
wide class of models, including Markov chains, Langevin 
equations, and deterministic systems connected to heat 
baths [7-10]. 

When one wants to test the validity of eq. (1) for a given 
system, one has to know AG12. From experimental point 
of view, this is sometimes difficult. Instead, eq. (1) can 
be indirectly verified by changing the magnitude of the 
external force / applied to the system. If the magnitude 
of the external force is changed by a small amount 5f, 
AG12 would be modified as AG12 AG12 -\-Sf where 
£ is a vector connecting mi and m2. Then, according to 
eq. (1), the ratio of the kinetic rates would be modified as 



k2l{f + Sf) k2l{f) 



exp {(35 f ■ I) . 



(3) 



ki2{f + df) kuif) 

Note that this condition can be verified without determin 
ing AG12. 



T. Harada 1 N. Nakagawa 2 



In the case of a molecular motor, the motor exhibits 
stepwize motion along a one-dimensional filament with a 
fixed step size It would be possible to view this be- 
havior as transitions among discrete states. Then, one 
would expect the condition given in eq. (3) holds for the 
kinetic rates of the forward and backward steps of the 
motor. Actually, it has been used in one-dimensional hop- 
ping models as a constraint for the transition rates in the 
model [11-13]. However, single molecule-experiments re- 
vealed that eq. (3) does not hold for a variety of molecular 
motors [1-5]. 

For instance, Nishiyama et al. measured the kinetic 
rates of forward and backward steps of conventional ki- 
nesin, denoted by fc+ and fc_, respectively, for varying 
magnitude of the external force / applied to the motor in 
parallel to the filament. They obtained data implying 



k.{f + Sf) k.{f) 



(4) 



where is a dimcnsionlcss parameter, whose value was ~ 
0.4. This result has been confirmed by subsequent studies 
[2,3]. Furthermore, the value of 6 has been measured 
for several other species of molecular motors, including 
22S dynein (taken from Tetrahymena cilia) [4] and myosin 
V [5] , and the results in those cases also indicate that the 
value of is significantly less than unity. 

Conventionally, the dimensionless quantity 9 introduced 
above has been referred to as the ratio of the sum of 
the "characteristic distances" , d, of forward and backward 
steps to the actual step size i: d = 9£ [1,3, 16]. Thus, 
the above experimental results imply that the character- 
istic distance is less than the actual step size of molecular 
motors. However, the physical significance of the charac- 
teristic distance has not been clarified yet. 

These experimental findings imply that cq. (3), and 
thus eq. (1), cannot be applied to the kinetic rates of the 
forward and backward steps of molecular motors. This 
empirical observation, however, may not suggest that the 
LDB condition [eq. (2)], from which eq. (1) is derived, is 
not applicable to microscopic models of a molecular mo- 
tor. Actually, through investigations of the Brownian mo- 
tor models [14,15], it has been realized that the kinetic 
rates k± does not satisfy eq. (1). However, the mecha- 
nism of violation of cq. (1) in these models has not been 
fully elucidated. This is because it is in general not easy 
to obtain simple expressions of the kinetic rates, k± , for 
the conventional Brownian motor models. Furthermore, 
in the earlier studies of the Brownian motor models, the 
LDB condition has not been always implemented. 

In this letter, we present a simple model of a molecu- 
lar motor that satisfies the LDB condition and for which 
k± can be explicitly calculated. For this model, we 
demonstrate that the dimcnsionlcss parameter 9 intro- 
duced above actually deviates from unity, and we clarify 
the mechanism of this deviation. This also clarifies the 
physical significance of the characteristic distance of a 



molecular motor. Finally, we demonstrate that the dimen- 
sionless parameter 9 provides a good measure to indicate 
the structure of the paths in the state space of the system. 

We here study a model possessing a "hidden" degree of 
freedom that is neither observed nor controlled experimen- 
tally. For example, this hidden variable could represent 
the extent of ATP hydrolysis reaction, i.e., the degree of 
completion of the reaction, which is not monitored in most 
experiments. This "hidden" degree of freedom is also con- 
ceptually equal to the internal state that is often incorpo- 
rated into a Brownian motor model. Here, for simplicity, 
we assume that there is only one such degree of freedom, 
denoted by y, and that its timescale is comparable to that 
of the spatial degree of freedom, x. 

Employing the above-stated assumptions, we study the 
following set of Langevin equations: 

= _~(!))^ ; + ,5) 
rur> = _^»W)+f + H(,). (6) 

Here, 7 and F are the friction coefficients for x and y, ^{t) 
and are zero-mean white Gaussian noises whose vari- 
ances arc 2"fk-QT and 2Tk^T, respectively, / represents an 
external force applied to the particle, and F represents a 
driving force conjugate to y. We impose a certain transla- 
tional symmetry on the potential U{x,y), which is speci- 
fied below, so that the motion of the particle is ^-periodic 
when projected onto the x axis. We also assume that the 
motion of the particle is L-periodic when projected onto 
the y axis, noting that the ATP hydrolysis reaction occurs 
in a cyclic manner, and we impose an appropriate trans- 
lational symmetry on the potential. With this interpreta- 
tion, the product of the driving force, F, and the periodic- 
ity of the system along the y direction, L, is interpreted as 
the change in the chemical potential occurring in a single 
ATP hydrolysis reaction: A/x = FL. Then, the total free 
energy is defined as V{x, y) = U{x, y) — fx—Fy. This type 
of model has been investigated as a model of molecular 
motors in several contexts [17,18]. Furthermore, by inter- 
preting y as an internal state of the motor, this model can 
be connected to standard Brownian motor models with 
discrete internal states [14, 15] by considering appropriate 
form of the potential U{x,y) and letting r/7 ^ 1. An im- 
portant property of the system described by eqs. (5) and 
(6) is that it satisfies the LDB condition [eq. (2)]. This 
can be easily confirmed by use of Onsager-Machlup path 
integral representation [9,19]. 

In order to grasp the general features of this model, we 
study several typical examples of the potential U{x,y). 
We first consider a simple potential function: U{x, y) = 
— UoCOs{TTx/£)cos{TTy/L) (see fig. la), which we term 
Model I. This potential possesses the reflection symme- 
try with respect to x: U{x,y) = U{xa — x,y) with 
^xq G [0,£). Let rm = (a;m,ym) denote a minimum of 
the energy function V{x,y). If the temperature is suffi- 
ciently low relative to the height of the barriers, AV, the 
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Fig. 1: a) Form of tlie potential U(x,y) = — (7o cos(7ra::/£) cos(7rj//L) (Model I) represented by a grayscale, with the value of U 
decreasing as the shade becomes darker. The small white dot labeled m denotes a minimum, and the small black dots labeled 
A, B, C and D denote saddles. The dashed lines (solid lines) associated with the fixed points indicate the stable (unstable) 
directions with positive (negative) eigenvalues of the Hessian matrix, Ti. b) The reversibility parameter as a function of the 
driving force F with / = 0. The solid curve represents the value of 9 calculated using eq. (12). The solid triangles, squares and 
circles represent the values of 9 calculated through numerical simulation of the Langevin equations (5) and (6) for /3f/o — 5, 10 
and 20, respectively. The error bars are smaller than the symbols, c) Schematic representation of the shift of the saddles under 
a perturbation, SE — {0,5F) (indicated by small arrows). The arrows SrA, Svb, Src, and 5ro represent the directions in which 
the saddles are shifted, as calculated with eq.(16). 



particle will hop from this minimum to one of its neigh- 
boring minima, + {£,L), rm + {£, -L), rm + {-i,L) 
or Tin + (— £, — L), through the saddle ta — (xajVa), 
rB = (a;B,2/B), t"c = {xc,yc) or ru = (xd.Z/d), respec- 
tively (see fig. la). Because the two minima + {£, ±L) 
located in the forward direction possess the same value of 
X, hops to them cannot be distinguished when we consider 
the particle motion projected onto the x axis. The same 
is true for hops in the backward direction. Note that the 
hopping motion of the particle becomes well described as 
a Markovian process in the low-temperature limit, even if 
one degree of freedom is not monitored. This is because 
whether the particle hops in the forward or backward di- 
rection becomes independent of all previous steps in the 
low-temperature limit. 

Using multi-dimensional Kramers rate theory [20-22], 
we find that the hopping rates from the well around to 
the wells around r,j-, + {£, L), -\- {£, —L), + {—£, L) 



and 



— -L) are given by 



«(+,+) 



Ab 



-^A/|detHm7^B'l e-'^i^^'"^^-''^'"^ (8) 
zvr 



k 



(-,-) 



'«( + ,-)C , 

K(+_+)e 



(9) 
(10) 



Here, Hm, T^a and Hb denote the Hessian of V{x, y) evalu- 
ated at rjii, r"A and Tb, respectively. The values Aa and Ab 
are the positive eigenvalues of —aTi.A and —ctHb, where a 
is a matrix whose components are an — I/7, (T22 — V-T, 
and ai2 = (J21 = 0. The periodicity of the system implies 
the equivalence of the saddles ta and ru and of the sad- 
dles rB and rc and that we have the relations a:D ~ xa—£, 



J/D = J/A - L, xc = XB - £, and yc = i/b + L. Then, eqs. 
(9) and (10) follow. 

The kinetic rates for forward and backward steps are 
given by k± = fc(±,+) + fc(±,-)- Therefore, from the defi- 
nition of 6 given in eq. (3), it can be expressed as 

_ XAk{+,+) + a;Bfc(+,-) a:cfc(_ + XY)k(_ _) 



+) 



(11) 



Thus, it is found that the value of 9 is given by an average 
of the X coordinates of the saddles in which the weights 
are the hopping rates. This is one of the main results of 
this letter. 

Because we are interested in the low-temperature limit, 
the hopping rates in the forward direction, and 
differ greatly in general. Let x-). denote the x co- 
ordinate of one of the two saddles in the forward direction 
at which the energy is lower. We term such a saddle a 
"dominant saddle". The situation is similar for the hop- 
ping rates in the backward direction, A;(-,-f) and fc(_.„). 
Let X- denote the x coordinate of the dominant saddle in 
the backward direction. Then, in the limit 13/S.V — > 00, we 
have 

9 = {i+-5:-)/£ + 0{T//W). (12) 

From this expression, it is found that the value of does 
not depend on the temperature if it is sufRciently low. This 
observation is consistent with recent experimental data [3] . 
It should be noted that in the case the dominant saddles 
are not uniquely determined (e.g., in the case f — F = 
for Model I), eq. (12) is not apphcable, while eq. (11) 
remains valid. 
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Using eqs. (11) and (12), we can calculate the value of 
6 once the potential U (x, y) is given. Figure lb displays 
the result of our numerical calculation for Model I with 
/ = 0. It is observed that the value of is equal to unity 
when the driving force, F, is vanishing. Moreover, it is 
found that 9 is always unity if F = 0, not only for the 
case / = but also / 7^ (data not shown; see the discus- 
sion below). However, if F is switched on, the value of 6 
deviates from unity. In Model I, we find that 9 decreases 
monotonically with the magnitude of F. We note that 6 is 
an even function of F because of the reflection symmetry 
of the system. In fig. lb, as a function of F, calculated 
through direct simulation of eqs. (5) and (6) is plotted for 
various values of the temperature. It is seen that the form 
of 9 seems to display the limiting behavior expressed in 
eq. (12) in the T ^ limit. 

We now present arguments that account for the obser- 
vations discussed above. First, the reason that we have 
9=1 when F = can be understood as follows. The 
transition probability of a microscopic path satisfies the 
LDB condition given in eq. (2), where G{r) = V{x,y). 
From this condition, by integrating out paths connecting 
ri = {xo,yo) and r2 = {x,y), one can prove that the tran- 
sition probability P{x,y,t\xo,yo,to) satisfies the relation 



P{x,y,t\xo,yo,to) 
P{xo,yo,t\x,y,tQ) 



= Q0[Vixo,Vo)-V{x,y)]_ 



(13) 



Because we observe the motion of the particle projected 
onto the x axis, the relevant transition probability is the 
projected probability 



P{x,t\xo,0) 



J dyo J dyPinit{xo,yo)P{x, y, t\xo, yp, 0) 
-Pinit(a;o) 



. (14) 

where Finit(a;,y) is the initial distribution and Pinit(a^) = 
J dyPinitix^y) represents a projected initial distribution. 
Because the probability distribution relaxes inside a po- 
tential well very quickly before the occurrence of a hop, 
the projected transition probability is insensitive to the 
choice of the initial distribution if t is much larger than 
the intra- well relaxation time, Tin. Considering this point, 
we choose Pinit{x,y) = C exp[—pU{x,y)], where C is a 
normalization constant. 

Then, using eqs. (13) and (14), for F = 0, we obtain the 
expression 



P{xm + e, t\x„„ 0)/P(a;m, t\xn, +£,0)= e' 



fife 



(15) 



where the identity Pinit(2;m) = Pinit{x^+£) has been used. 
Next, these transition probabilities are related to the ki- 
netic rates as P{xm ± i,At\xm,0) = k±At/fl + 0{At^) 
for At satisfying <C At <C {k±)~^, where fi = 
/ dxdy exp[—(3V{x, y)] and A denotes the interior of the 
potential well around rm [23]. Then, eq. (15) results in 
k+/k- = ex.p{(3fi), and ^ = 1 is obtained for any value of 

Next, we explain why ^ is a decreasing function of |F|. 
The positions of the saddles, from which the value of 6 is 



estimated according to eq. (12), are solutions of the fixed 
point equations dV{x,y)/dx = and dV{x,y)/dy = 0. 
Let rx(-E) = (xxiE) , yx{E)) , where E = (/,F), denote 
such a solution. When the force E is slightly modified as 
E E + SE, where SE = {Sf, SF), the fixed point moves 
as 

rxiE) ^ rxiE+SE) = rx{E)+n^^SE+OiSE^), (16) 

where Hx is the Hessian of V{x,y) evaluated at rx{E). 

The shifts undergone by the saddle points due to a 
change in E, as given by eq. (16), arc depicted in fig. Ic. 
In the case of a positive driving force (F > 0), the domi- 
nant saddles are ta and rc- In Model I, the eigenvector 
of TYa at ta with a negative eigenvalue always lies in the 
first or third quadrant. For this reason, eq. (16) leads to 
the inequality xaIE + SE) - xaIE) < 0{SE'^) for 6f = 
and SF > 0. Analogously, we find xcAE + SE) - xc{E) > 
0{SE^). A slight increase of F thus results in a decrease 
of e, in other words, we have de{SE)/dF < for F > 0. 
Therefore, for F > and fixed /, 9 is a. monotonically 
decreasing function of F. A similar argument can be ap- 
plied in the case of negative F, and we find that for F < 
and fixed /, is a monotonically decreasing function of 
|F|. In this way, the shift of the dominant saddles and, 
thus, the change in the value of 6 are determined by the 
configurations of the eigenvectors at the saddles. 

Summarizing the above arguments, in the case of Model 
I, we can conclude that 9 decreases from unity as the mag- 
nitude of F increases from zero. 

In order to check the generality of the above argument, 
we have investigated several other models. Here, we con- 
sider one without the reflection symmetry with respect to 
x: U{x,y) 7^ U{xo — x,y) for '^xq S [0,^). The potential 
profile is plotted in fig. 2a (to which we refer as Model 
II) . Because of the lack of the reflection symmetry in this 
model, the particle will realize a non-zero mean velocity in 
the X direction even in the absence of the external force /, 
provided that the driving force F is non-zero. This situa- 
tion is reminiscent of that for a biological molecular motor 
(see the upper inset of fig. 2b). 

In fig 2b, 9 calculated using eq. (12) is plotted as a func- 
tion of F with / = 0. It is found that 9 = 1 over a certain 
range of values of F. Investigation of the potential at the 
saddles reveals that within this range of F, the dominant 
saddles are ta and Td- Therefore, by considering eq. (12) 
and the identity xu = xa — we conclude that 9 = 1. 
Note that in this case, the dominant path for a forward 
step is the reverse of that for a backward step. Conse- 
quently, the trace of many sequential steps remains in a 
one-dimensional sub-space (see inset 1 of fig. 2b). Let us 
term this situation "dominant-path reversible" . 

However, at a certain critical value of F, denoted by 
Fc, 9 drops discontimiously from unity, and for F > F^, 
is a continuously decreasing function of F. At F = Fc, 
the roles of saddles are exchanged; the dominant saddle 
in the backward direction changes from ru to rc- Thus, 
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Fig. 2: a) Form of the potential i7(a;, j/) = —Uo[l + acosn {x/£ — y /L)][l + cos (nx/l) cos {ny / L)], with a = 0.5 (Model II), 
represented by a grayscale, with the notation the same as in fig. 1. b) The reversibility parameter S as a function of F with 
/ = 0. The solid curve represents the value of calculated using eq. (12). The structures of the dominant paths in the state 
space are depicted in insets 1 and 2 for each branch of the curve. The upper inset displays the mean velocity as a function of 
the driving force F with f = for (3Uo = 100. The vertical axis is in units of 10'^'^ ksT /{-(£). c) 9 as a function of / with 
F — Uo/L calculated using eq. (12). 



for F > Fc, the dominant saddles are ta and rc, which 
become closer as F increases for the same reason as that 
for Model I. Then, from eq. (12), it is concluded that the 
value of 9 decreases as F increases. In this case, the dom- 
inant path for a backward step is not the reverse of that 
for a forward step, and the trace of sequential steps forms 
a branched structure in the state space (see inset 2 of 
fig. 2b). In contrast to the dominant-path-reversible situ- 
ation, we term this situation "dominant-path irreversible" . 
Note that a similar phenomenon occurs for negative F be- 
cause of the symmetry of the system. Moreover, this type 
of switching also occurs when the load, /, is varied with 
fixed F, as plotted in fig. 2d. 

In this way, the value of 6 directly reflects the structure 
of the dominant paths in the state space. It is interesting 
that the conventionally analyzed quantities such as the 
mean velocity and the diffusion cocfflcient do not exhibit 
such a sharp transition. Thus, it is the value of 0, not 
these other quantities, that provides the information of 
the structure of the dominant paths in the state space of 
the system. 

In conclusion, the physical significance of the dimen- 
sionless parameter 9, which has been defined as the ratio 
of the characteristic distance to the actual step size of 
a molecular motor, has been investigated by considering 
simple Langevin models that have a hidden degree of free- 
dom. We have identified the mechanism through which 9 
deviates from unity in these models. It was found that if 
no driving force is applied to the hidden degree of free- 
dom, 9 is always unity. This means that if we find 6 1 
in a certain system, there could be at least one hidden 
degree of freedom that is directly coupled to an external 
energy source. Furthermore, 9 is always unity even when a 
driving force is present, if the system is in "dominant-path 
reversible" situation, i.e., the dominant path for a forward 
step is the time-reverse of the dominant path for a back- 
ward step. By contrast, 9 deviates from unity if the system 



is in "dominant-path irreversible" situation, i.e., the dom- 
inant path for a forward step is not the time-reverse of the 
dominant path for a backward step. In this case, the dom- 
inant paths have a branched structure in the state space. 
Thus, the value of 9 serves as an indicator of the structure 
of the dominant paths in the state space even when several 
degrees of freedom are not monitored. These features of 9 
are expected to be valid if the number of hidden degree of 
freedom exceeds one, since the main results presented in 
this letter are easily extended to systems with more degree 
of freedom. 

On the basis of the above results, we now present an 
interpretation of relevant experimental findings for molec- 
ular motors with regard to dominant-path reversibility. 
As noted above, it has been found that most of biological 
motors exhibit values of 9 less than unity. According to 
our discussion, there could be hidden degrees of freedom 
that are coupled to the energy source (i.e. the reservoir 
of ATP, ADP and Pi) and that the dominant paths in 
the state space are branched. This seems to imply that 
these motors in general hydrolyze ATP molecules even in 
backward steps (dominant-path irreversible) . This is con- 
sistent with several experimental results on conventional 
kinesin [1,2]. 

In Model II, there was a situation where 9 remains unity 
while non-zero velocity is exhibited (dominant-path re- 
versible). In this case, the dominant path for a backward 
step is the reverse of the dominant path for a forward 
step, which implies that ATP is synthesized when the mo- 
tor steps backward. This situation could correspond to 
the case of Fi ATP synthase that synthesizes ATP when 
it is forced to rotate backward [24,25]. Therefore, it is 
expected from our argument that the value of 9 for this 
molecule is close to unity. Although the value of 9 has not 
been determined for this molecule, its determination will 
serve as an experimental test of the present theory. 

To properly interpret the results of this work, it should 
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be noted that the connection between the models studied 
here and actual molecular motors is not clear, as these 
models are quite simple. For instance, the actual molecu- 
lar motors are expected to have several substates those are 
not necessarily coupled to its stepwise motion. Since such 
substates can also be incorporated to the model studied in 
this letter by modifying the profile of the potential U{x,y), 
one would be able to see the effects of such substates on 
the property of the characteristic distance. Moreover, in 
most models we have studied, including two models exem- 
plified here, the value of 6 is less than or equal to unity. It 
is not clear at present whether there is a case in which 9 is 
greater than unity. Further theoretical and experimental 
studies of the parameter 6 should lead to a better under- 
standing of the mechanisms and architectures of biological 
molecular motors. 
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